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Abstract 

Post-transcriptional regulatory mechanisms are of fundamental importance to form robust genetic networks, but their roles 
in stem cell pluripotency remain poorly understood. Here, we use freshwater planarians as a model system to investigate 
this and uncover a role for CCR4-N0T mediated deadenylation of mRNAs in stem cell differentiation. Planarian adult stem 
cells, the so-called neoblasts, drive the almost unlimited regenerative capabilities of planarians and allow their ongoing 
homeostatic tissue turnover. While many genes have been demonstrated to be required for these processes, currently 
almost no mechanistic insight is available into their regulation. We show that knockdown of planarian Notl, the CCR4-N0T 
deadenylating complex scaffolding subunit, abrogates regeneration and normal homeostasis. This abrogation is primarily 
due to severe impairment of their differentiation potential. We describe a stem cell specific increase in the mRNA levels of 
key neoblast genes after Smed-notl knock down, consistent with a role of the CCR4-N0T complex in degradation of 
neoblast mRNAs upon the onset of differentiation. We also observe a stem cell specific increase in the frequency of longer 
poly(A) tails in these same mRNAs, showing that stem cells after Smed-notl knock down fail to differentiate as they 
accumulate populations of transcripts with longer poly(A) tails. As other transcripts are unaffected our data hint at a 
targeted regulation of these key stem cell mRNAs by post-transcriptional regulators such as RNA-binding proteins or 
microRNAs. Together, our results show that the CCR4-N0T complex is crucial for stem cell differentiation and controls stem 
cell-specific degradation of mRNAs, thus providing clear mechanistic insight into this aspect of neoblast biology. 
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Introduction 

Post-transcriptional control is central for the regulation of gene 
expression in stem cells [1]. A key post-transcriptional process is 
mRNA degradation [2-4] the control of which is believed to be as 
important as transcriptional regulation [5,6]. Although transcrip- 
tional regulation has been extensively studied, less is known about 
the developmental and physiological roles of mRNA degradation 
in stem cells, which are thought to involve the same RNA binding 
proteins [7] that act together to coordinate many complex aspects 
of mRNA biology, one of which is degradation. 

mRNA degradation starts with deadenylation (i.e. shortening of 
the poly(A) tail) [2,8]. This affects gene expression both by 
decreasing translational activity and mRNA stability [9,10]. The 
major deadenylase in eukaryotes is the CCR4-NOT complex [1 1— 
13], which is also involved in regulating several other aspects of 
mRNA metabolism, such as mRNA export [14,15], translation 
[15] and transcription itself [11,16-18]. 



In yeast, the CCR4-NOT complex is composed of nine different 
subunits [11]: two deadenylases (Ccr4p and Pop2p/Caflp), five 
Not proteins (Notlp-Not5p), Caf40p and Cafl30p. Among them, 
Notlp, a 240 kDa protein, is thought to act as a scaffold and is the 
only subunit required for yeast viability [11,19]. Most of the 
subunits of the yeast complex are conserved across metazoans [20- 
22], with the exception of Not5p and Cafl30p. In mammals two 
paralogous genes with mutually exclusive expression patterns 
encode each deadenylase of the complex [23]. Furthermore, the 
two deadenylase subunits Cafl and Ccr4 regulate distinct sets of 
mRNAs [24,25]. 

A number of translational repressors interact with the CCR4- 
NOT complex to repress their targets. For instance, Nanos 
proteins [26,27], PUF proteins [28,29], Smaug [30], and Bicaudal- 
C [3 1] all repress their target mRNAs via interaction with different 
subunits of the CCR4-NOT complex. Furthermore, the CCR4- 
NOT complex mediates the deadenylation of miRNA-targeted 
and piRNA-targeted mRNAs, executing the repressive functions 
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Author Summary 

Although transcriptional regulation in stem cells is a very 
active subject, much less is known about how post- 
transcriptional mechanisms of gene expression affect stem 
cells. Here, we use freshwater planarians in order to 
address this question. Planarians have a striking regener- 
ative capacity driven by a population of pluripotent stem 
cells, the neoblasts. Control of both proliferation and 
differentiation is thought to rely heavily on post-transcrip- 
tional mechanisms, but their precise role is unknown. 
Poly(A) tail length regulation is an important mechanism of 
post-transcriptional control of gene expression as changes 
can be very rapid, and longer poly(A) tails are linked to 
increased mRNA stability and translational activity. We 
investigated the role of the CCR4-N0T complex, the major 
deadenylating complex in eukaryotes, by knocking down 
its main scaffolding subunit called Notl. Neoblasts in 
knock down animals are unable to differentiate and 
accumulate mRNAs with longer poly(A) tails. Our results 
show that the CCR4-NOT complex is needed for the 
targeted degradation of specific mRNAs expressed in stem 
cells, and the failure of this process likely prevents 
neoblasts from differentiating. These results reveal a new 
functional aspect of the CCR4-NOT complex and offer a 
mechanistic insight into the regulation of planarian stem 
cells. 

of some small RNAs [32-35]. The CCR4-NOT complex directly 
binds to GW182, a component of miRNA repression complexes 
through evolutionary conserved motifs [36,37]. 

Little is known, however, about the role of the CCR4-NOT 
complex in stem cells. It was found, for instance, that different 
components are important in maintaining mouse and human ESC 
identity [38] , but the mechanisms remain largely unexplored. The 
freshwater planarian Schmidtea mediterranea is an emerging model for 
stem cell biology [39-41]. Its striking regeneration capacities are 
sustained by the presence of the neoblasts, a population of 
pluripotent stem cells that not only drive regeneration but sustain 
constant homeostatic cell turnover as well [42]. Planarian 
neoblasts can be eliminated by irradiation and are amenable to 
RNAi-mediated gene knock down. Furthermore, their abundance 
in the organism allows the quantitative evaluation of phenotypes. 

The regulation of neoblasts and their pluripotency involves both 
transcriptional and post-transcriptional regulation, and a number 
of putative post-transcriptional regulators have been described to 
affect neoblast function [43-46]. Neoblasts also contain chroma- 
toid bodies, RNA rich granules similar to germ granules which are 
thought to constitute hubs for mRNA processing and post- 
transcriptional regulation [39,45,47-49] . Recently, transcriptomic 
profiles of neoblasts have become available [50-54], confirming 
their long known resemblance to germ line stem cells [55], but also 
highlighting the importance of post-transcriptional mechanisms 
for their regulation [50] and the conservation of pluripotency 
determinants between planarian and mammalian stem cells 
[51,52,56,57]. 

Here, we use S. mediterranea to investigate the role of the CCR4- 
NOT complex in stem cell regulation by characterizing the 
function of the Smed-notl gene, which encodes for the homologue 
of Notlp in yeast and CNotl in mammals. We report that its 
knock down specifically affects the differentiation and self- 
maintenance capabilities of neoblasts. Smed-notl knock down 
results in a progressive increase in the levels of several neoblast 
transcripts, and we demonstrate that this increase is stem cell 
specific. Finally, we observe that these same mRNAs have stem 



cell specific increases in the frequency of long poly(A) tails after 
Smed-notl knock down, showing that the observed increases in 
mRNA levels in stem cells are likely a consequence of decreased 
targeted deadenylation by the CCR4-NOT complex. Our findings 
highlight a likely central role for poly(A) tail length regulation in 
orchestrating pluripotent stem cell differentiation. 

Results 

In silico characterization of the CCR4-NOT complex in S. 
mediterranea 

We identified the different subunits of the CCR4-NOT complex 
in the planarian species S. mediterranea by TBLASTN searches in 
the S. mediterranea genome and in our reference transcriptome 
assembly [50,53] and other genomic and transcriptomic resources 
[58-60] (Table SI). We identified homologues of all metazoan 
genes known to encode subunits of the CCR4-NOT complex. 
Two components of the CCR4-NOT complex, the deadenylases 
Ccr4 and Cafl were previously described in the planarian species 
Dugesia japonica [45] . However, no phenotype was reported for 
these enzymatic components after RNAi-mediated knock down 
and we also observed no strong phenotype (Figure SI) for Smed- 
not6 (the Ccr4 orthologue, Table SI), Smed-not7A and Smed-not7B 
(the two paralogues of Cafl in S. mediterranea, Table SI). We 
instead chose to focus on Smed-notl, as Notl is believed to act as 
the central scaffolding protein in the complex and it is the only 
component of the complex essential for viability in yeast [11,19]. 

Smed-notl is expressed in CNS and throughout the 
parenchyma in an irradiation-sensitive manner 

We investigated the expression pattern of Smed-notl by whole 
mount in situ hybridization (WMISH). We observed broad 
expression of this key component of the core deadenylation 
complex, including expression throughout the parenchyma and 
the central nervous system (CNS) (Figure 1A, top panel). This 
pattern suggested to us that Smed-notl may be expressed in 
neoblasts, since they are distributed throughout the parenchyma. 
To test this we monitored parenchymal expression after irradiation 
to remove neoblasts and observed that the parenchymal compo- 
nent of expression disappeared over a period of 5 days after 
irradiation (Figure 1A, middle and bottom panels). AH neoblasts 
disappear 24—48 hours after lethal irradiation, and consequendy 
the expression of neoblast specific genes disappears over a similar 
period. As controls we analyzed the expression of Smedtud-1, the 
orthologue of the previously described Schmidtea polychroa Tudor 
gene Spoltud-1 [48] (Figure IB), Smed-va.sa-1, a Vasa orthologue of 
S. mediterranea [54,61] (Figure 1C), and Smed-pcna, the orthologue of 
the PCNA gene described in Dugesia japonica [62] (Figure ID). As 
expected, the neoblast-specific signals of all three disappeared by 
day 3 post-irradiation, while irradiation insensitive expression in 
differentiated cells of the CNS remained for Smedtud-1 and Smed- 
vasa-1 (Figure 1B-D). Smed-eye53, a marker expressed in differen- 
tiated cells [63], was used as a control to demonstrate that gene 
expression in post-mitotic cells is not ablated by irradiation 
(Figure IE). Since Smed-notl hybridization signals present in the 
parenchyma were reduced but did not completely disappear by 
day 3 post-irradiation, it is likely that Smed-notl is expressed in 
neoblasts and their recent progeny, but also other post-mitotic 
cells. These data indicates that Smed-notl is broadly expressed 
throughout the planarian body, consistent with a housekeeping 
function, in a pattern that includes neoblasts, their progeny as well 
as differentiated cells. Investigation of Notl expression in previous 
transcriptome based studies of mRNAs expressed in neoblasts is 
consistent with this expression pattern (Figure S2). 
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Figure 1. Smed-notl is expressed in CNS and throughout the parenchyma in an irradiation-sensitive manner. (A-D) WMISH of Smed- 
notl (A), Smedtud-1 (B), Smed-vasa-1 (C), Smed-pcna (D) and the control marker Smed-eyeS3 (E), in non-irradiated and irradiated animals. Smed-notl 
signals are detected as a broad staining pattern in non-irradiated animals (A, top panel, non irrad), and decrease 3 days (A, middle panel, irrad 3d) and 
5 days (A, bottom panel, irrad 5d) after lethal irradiation. Smed-notl signal is detected in the CNS (arrows). The neoblast specific signals of other 
mRNAs expressed in neoblasts disappear by day 3 of irradiation (B, C, D, irrad 3d), while signals are still detectable in the CNS for Smedtud-1 (B, 
arrows), and Smed-vasa-1 (C, arrows). No Smed-pcna signal is detected in the CNS (D). The differentiated cell control marker Smed-eye53 shows no 
differences upon irradiation (E). Anterior is to the left. Scale bars: 500 |j.m. 
doi:1 0.1 371 /journal.pgen.1 004003.g001 



Smed-notl is required for planarian regeneration and 
homeostatic cell turnover 

We then analyzed the function of Smed-notl by RNAi 
experiments. All Smed-notl (RNAi) animals displayed abrogated 
regeneration capacities and eventually died. They were able to 
produce both anterior and posterior regeneration blastemas, but 
never completed the regenerative process (Figure 2B, vs. 
Figure 2 A). In order to test if the formation of a regeneration 
blastema depended on the time of transection, we cut animals at 1 , 
3, 5, 10 and 15 days after Smed-notl (RNAi) treatment (Figure S3). 
We found that animals were able to produce a regeneration 
blastema at all-time points, however animals cut earlier produced 
larger blastemas. All blastemas of Smed-notl (RNAi) worms eventu- 
ally regressed. The ability of Smed-notl (RNAi) animals to produce a 
large regeneration blastema at early time points after RNAi 
suggests that mitotic neoblasts, the source of blastema cells, are still 
present and proliferating. 

We then analyzed the phenotype of intact Smed-notl (RNAi) 
animals. We found that Smed-notl (RNAi) worms presented head 
regression (Figure 2D, vs. Figure 2C) [43,44,64,65], a symptom of 
interrupted homeostatic cell turnover. A survival curve during 
which the onset of tissue homeostasis defects was recorded (N = 40) 
demonstrated temporal phenotypic variability (Figure 2E). In 
Smed-notl (RNAi) animals these defects were seen first at 15 days 
after RNAi, and in the majority of animals after 20 days 
(Figure 2C). Variable degrees of head regression were also 
observed after 20 days of RNAi (Figure 2C, also see Figure S3). 
By day 22 all animals had defects, showing complete penetrance of 



Smed-notl RNAi. All animals died by day 36 after dsRNA delivery, 
with the majority of deaths occurring between day 26 and day 34 
(Figure 2E). All control(RNAi) animals survived without any defect 
for >35 days, however. Together, these results demonstrate that 
Smed-notl is needed for regeneration and homeostatic cell turnover 
in S. mediterranea. 

Smed-notl(RNAi) animals maintain proliferative neoblasts 

Next, we analyzed the mitotic marker phospho-histone-H3 
(h3p) in Smed-notl (RNAi) animals at 5, 10, 15 and 20 days 
after RNAi. Up to 15 days all (N = 7 per time point) had 
mitotic neoblasts comparable in numbers to those of 
control(RNAi) animals (Figure 3A). After 20 days all animals 
still had mitotic cells, but at variable levels, in agreement with 
our earlier phenotypic characterisation (Figure 3B-E). Most 
had normal levels (Figure 3C, vs. Figure 3B), although 
animals with more severe head regression defects showed a 
visible reduction in the mitotic levels (Figure 3D-E, vs. 
Figure 3B-C), but overall the reduction in mitoses was not 
statistically significant. These experiments show that Smed- 
notl (RNAi) worms have mitotic cells, even as head regression 
defects progress. Significantly, similar defects are seen in 
irradiated worms only weeks after complete loss of mitotic 
activity. Therefore, we interpret our data as showing that 
effects on stem cell proliferation are not the primary cause of 
the Smed-notl (RNAi) phenotype, instead implicating neoblast 
differentiation impairment as responsible for regenerative 
failure, head regression and other defects. 
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Figure 2. Smed-not1 is required for planarian regeneration and homeostatic cell turnover. (A-B) Control(RNAi) (A) and Smed-notl(RNAi) (B) 
animals cut 5 days after RNAi, and monitored every 2 days after transection. Top panels: anterior wounds; bottom panels: posterior wounds. Smed- 
notl animals are able to produce blastemas (white tissue) in both anterior and posterior wounds but fail to regenerate. (C-D) Intact control(RNAi) (C) 
and Smed-notl (RNAi) (D) animals 20 days after RNAi, anterior side. Smed-notl animals 20 days after RNAi display variable levels of head regression 
defects. (E) Survival curve of Smed-notl animals (N=40), indicating the percentage of animals without any defects (grey) and the percentage alive 
(blue). All control(RNAi) animals survived without defects for more than 40 days (N = 40). Scale bars: 500 jj.m. See also Figure S3. 
doi:1 0.1 371 /journal.pgen.1 004003.g002 



By performing FACS experiments we found that Smed- 
not(RNAi) led to a moderate reduction of the sorted irradiation 
sensitive XI cells, which primarily contains neoblasts, at 15 
days but not at 10 days (Figure 3F. Figure S4A-B). Even 
though no significant decrease of h3p cells was detected at this 



time point we interpret our FACS data as more sensitive and 
conclude that both methods consistently detect large numbers 
of proliferating neoblasts 15 days after Smed-notl dsRNA 
administration, further implicating neoblast differentiation 
defects instead. 
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Figure 3. Smed-notl(RNAi) animals maintain mitotic neoblasts. 

(A) Quantification of mitosis by counting of h3p-positive cells in whole 
mount immunohistochemistry on control(RNAi) and Smed-notl(RNAi) 
animals 5, 10, 15 and 20 days after RNAi (N = 7 per time point). No 
significant differences are detected. Representative control(RNAi) (B) 
and different Smed-notl(RNAi) worms (C-E) 20 days after RNAi, 
immunostained with the mitotic marker h3p (h3p, green) and 
counterstained with phospho-tyrosine (p-tyr, red) in order to show 
head regression defects. Smed-not1 (RNAi) animals still display detect- 
able mitotic cells, even as head regression defects occur. The number of 
mitotic cells detected is smaller in the animals with the most severe 
head regression defects (D-E). Anterior is to the left. Scale bars: 500 am. 
(F) Quantification of FACS sorted X1 cells in control(RNAi) and Smed- 
notl(RNAi) animals 10 and 15 days after RNAi, and wild type irradiated 
animals. While no significant differences are observed 10 days after 
knock down, Smed-notl(RNAi) animals show a reduced but significant 
decrease in percentage of XI cells. Error bars represent standard 
deviation and asterisks represent statistical significance. See also Figure 
S4. 

doi:10.1371/journal.pgen.1004003.g003 



Dynamics of neoblast cells and their progeny in Smed- 
notl(RNAi) animals 

In order to monitor the behaviour of neoblasts and their post- 
mitotic progeny during progression of the Smed-notl knock down 
phenotype we analysed the expression of neoblast and progeny 
markers [66]. Smedwi-1, a marker of neoblasts, Smed-nb.21.11e, a 
marker of early neoblast progeny, and Smed-agat- 1 , a marker of late 
neoblast progeny, were analyzed in control(RNAi) worms 
(Figure 4A-C) and Smed-notl (RNAi) worms after 10 (Figure 4D- 
F), 15 (Figure 4G-I), and 20 (Figure 4J-L, Figure S5A-F) days of 
RNAi. Only one time point (10 days) is shown for control(RNAi) 
worms, since no differences were observed between time points. 
Smedwi-1 expression was qualitatively the same after 10 and 15 
days of Smed-notl (RNAi), (Figure 4D, G, vs. Figure. 4A), but clearly 
reduced to a variable extent after 20 days (Figure 4J vs. Figure 4A, 
D and G; Figure S5A vs. Figure S5D); some animals had nearly 
normal expression while Smedwi-1 expression was severely reduced 
in those with the most severe head regression defects. These 
results, like those above, suggest that prominent stem cell loss 
occurs only when Smed-notl (RNAi) animals begin to regress the 
head and to die, again implicating differentiation impairment 
instead of proliferation or self-renewal as a primary cause for 
regenerative failure. 

Expression of Smed-nb.21.1 le, a marker of early neoblast 
progeny, looked broadly equivalent to that of control worms in 
animals fixed 10 days after Smed-notl dsRNA delivery (Figure 4E, 
vs. Figure 4B) but a clearly visible progressive decrease in the 
number of Smed-nb.21.1 le-positive cells was detected after 15 
(Figure 4H, vs. Figure 4B) or 20 (Figure 4K, vs. Figure 4B) days 
after RNAi. After 20 days of Smed-notl knock down animals had 
only a few remaining Smed-nb.21.1 le-positive cells (Figure S5E). 
These results show that clearly visible decreases in early progeny 
cell number precede the prominent decrease in neoblasts 
themselves. When we checked the expression of Smed-agat- 1 , a 
marker of late neoblast progeny, we observed a consistent 
qualitative increase in Smed-agat- 1 -positive cells in Smed-notl (RNAi) 
worms after 10 days (Figure 4F vs. Figure 4C). At later time points, 
however, the number of Smed-agat- 1 -positive cells also progressively 
declined (Figure 41, L, vs. Figure 4C). Again, we observed a 
considerable variability in Smed-notl (RNAi) worms after 20 days 
(Figure S5F). However, in all animals with a clear decrease in the 
number of Smed-agat- 1 -positive cells this was particularly apparent 
in the anterior region, a characteristic feature of Smed-agat-1- 
positive cell depletion upon neoblast elimination by irradiation or 
perturbation by RNAi [50,64,65,67] . 

Therefore, neoblasts are abundant 15 days after RNAi, and 
only clearly start to be depleted later, coinciding with the onset of 
head regression defects. Similar defects are observed after 
irradiation, however, these take >10 days to manifest after 
elimination of mitotic activity. In contrast Smed-notl (RNAi) animals 
display these defects when neoblasts are still present. These results 
suggest that a primary defect in neoblast differentiation, rather 
than neoblast maintenance, causes failure in tissue homeostasis. In 
support of this idea, altered stem cell progeny numbers precede the 
disappearance of Smedwi-1 signals and mitotic activity. This 
alteration can be observed as early as 10 days after RNAi in the 
case of Smed-agat- 1 -positive cells and 15 days for Smed-nb.21.1 le- 
positive cells, which are clearly depleted at this time point. 

To further test this, we compared the dynamics of the neoblast 
and progeny cell markers in Smed-notl (RNAi) animals to those of 
Smedwi-2 (RNAi) animals, in which neoblast differentiation is 
abrogated [44]. Smedwi-1, Smed-nb.21.1 le and Smed-agat-1 were 
expressed in Smedwi-2 (RNAi) animals with very similar dynamics to 
Smed-notl (RNAi) animals (Figure S5G-R). Taken together, these 
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Figure 4. Dynamics of neoblasts and their progeny in Smed-notl(RNAi) animals. (A-L) WMISH of the neoblast marker Smedwi-1 (A, D, G, J), 
the early neoblast progeny marker Smed-nb.21.1 1e (B, E, H, K) and the late neoblast progeny marker Smed-agat-1 (C, F, I, L) in control(RNAi) animals 
(A-C) and Smed-notl (RNAi) animals 10 (D-F), 15 (G-l) and 20 (J-L) days after RNAi. Smed-notl(RNAi) animals have detectable expression of Smedwi-1 
in all time points (A, D, G, J), although a decline in the level of Smedwi-1 signals is detected 20 days after RNAi (J). The dynamics of progeny markers is 
also abnormal, with a progressive decline of Smed-nb.21.1 1e signals (H, K) and an accumulation (F) followed by a decline (I, L) of Smed-agat-1 signals. 
Anterior is to the left. Scale bars: 500 um. See also Figure S5. 
doi:10.1371/journal.pgen.1004003.g004 



results show that, similar to Smedwi-2 RNAi, Smed-notl RNAi 
impairs neoblast differentiation with proliferation only affected in 
a later time point. 

Smed-notl (RNAi) animals have abnormal numbers of 
Smed-agat-1 transcripts and Smed-agat-1 -positive cells 

In order to achieve a quantitative measure of stem cell progeny 
mRNA levels in Smed-notl (RNAi) animals we performed quantita- 
tive real time PCR experiments (qRT-PCR) of Smed-nb.21.1 le and 
Smed-agat-1 transcripts in Smed-notl (RNAi) worms. We focused on 
earlier time points of 10 and 15 days after RNAi. qRT-PCR 
showed that Smed-nb.21.1 le levels in whole animals were similar to 
those of control(RNAi) after 10 and 15 days of dsRNA administra- 
tion (Figure 5A). However, Smed-agat-1 mRNA levels increased by 
day 10 and were almost two-fold higher after 15 days. Since this 
result did not correlate with what we observed by colorimetric 
WMISH, we quantified Smed-agat-1 -positive cells by fluorescent 
WMISH (FWMISH). A significant increase in the numbers of 
Smed-agat-1 -positive cells was found in Smed-notl (RNAi) animals 10 
days after dsRNA administration (Figure 5B), but their numbers 
declined to control levels after 15 days. These data confirmed the 
qualitative data from colorimetric WMISH (Figure 4). The 
distribution of Smed-agat-1 -positive cells in Smed-notl (RNAi) animals 
1 5 days after dsRNA delivery was different from controls, with less 
Smed-agat-1 -positive cells in the anterior part of the worm 
(Figure 5C). This revealed a stark discordance between the 
number of Smed-agat-1 -positive cells and the level of Smed-agat-1 
mRNA. We conclude that Smed-agat-1 transcripts are accumulat- 
ing in decreasing numbers of Smed-agat- 1 -positive cells after 15 
days of Smed-notl knock down, and that each Smed-agat-1 -positive 



cell contains an increased average number of Smed-agat-1 
transcripts. A similar, but less pronounced, process could also 
explain reduced numbers of Smed-nb.21.1 le positive cells and 
discordant stable levels of this transcript, which do not drop by day 
15. 

Given that the CCR4-NOT complex is known to regulate gene 
expression through its deadenylating activity we wished to 
ascertain if it could be directly responsible of the discordance 
between Smed-agat-1 positive cell number and mRNA levels. If an 
increase in mRNA levels is caused by impaired deadenylation and 
subsequent degradation we would expect to observe increased 
frequency of long poly(A) tail lengths. Using a poly(A) tail length 
(PAT) assay [68-71], in whole worms we observed that this was 
the case for both Smed-agat-1 and Smed-nb.21.1 l.e (Figure 5D), with 
the first giving starker differences, while control mRNAs Smed-eif-3, 
Smed-mhc, or Smed-ef-2 were only mildly affected or not affected at 
all. In addition a spike-in control of exogenous mRNA showed 
equal poly(A) tail length distribution across samples (Figure 5D), 
showing that the differences observed are present in our different 
planarian mRNA samples and not introduced by the PAT assay 
technique. 

Smed-notl(RNAi) animals have increased levels of 
transcripts expressed in stem cells with increased 
frequency of long poly(A) tails 

Next, we found that WMISH analysis of neoblast markers, 
Smedtud-1, Smed-vasa-1 and Smed-pcna, was suggestive of increased 
levels of these transcripts in Smed-notl (RNAi) worms, with 
qualitatively visible differences after both 10 and 15 days after 
RNAi treatment (Figure 6A— C). These data further demonstrate 
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Figure 5. Smed-notURNAi) animals have increasing numbers of Smed-agat- 7 transcripts with increased frequency of long poly(A) tails 
but decreasing numbers of Smed-agat- 7-positive cells. (A) Quantification of the level of expression by qRT-PCR of the progeny markers Smed- 
nb.21.l1e and Smed-agat-1 in Smed-notl(RNAi) animals 10 and 15 days after RNAi, normalized expression and relative to respective control(RNAi) 
samples. Error bars represent standard deviation, asterisks represent statistical significance. Smed-agat-1 transcripts accumulate progressively after 10 
and 15 days of RNAi. (B) Quantification of the number of Smed-agat- 7-positive cells in control(RNAi) and Smed-notl(RNAi) worms 10 and 15 days after 
RNAi. Animals (N = 8 per time point and treatment) were stained by FWMISH of Smed-agat-1 and analyzed by confocal microscopy. Multiple squares 
corresponding to 250 ^m 2 in both the dorsal and ventral parts of the animals were selected for counting along the length of the animals. 130 squares 
were counted, 4 representative squares are shown. Error bars represent standard error of the mean, asterisk represents statistical significance. A 
significant accumulation of Smed-agat- 1 -positive cells is detected in Smed-not1 (RNAi) animals 10 days after RNAi but a decrease is observed at 15 days 
(C) Confocal Z-projection of FWMISHs of Smed-agat1 in control(RNAi) (top panels) and Smed-notl(RNAi) (bottom panels) animals corresponding to 
anterior dorsal (left) and anterior ventral (right) regions of representative worms. Smed-agat- 1 -positive cells are shown in green, nuclei 
counterstaining in magenta. Scale bar: 100 um. (D) PAT assays reflecting the distribution of mRNA poly(A) tail lengths for Smed-agat-1 and Smed- 
nb.21.11e, the housekeeping mRNAs Smed-ef-2 and Smed-eif-3, the tissue specific mRNA Smed-mhc and a spiked-in control mRNAs in control(RNAi) (c) 
and Smed-not1 (RNAi) (n) animals 10 and 15 days after RNAi. Size markers used are represented in blue, the theoretical length of the amplicon 
corresponding to the deadenylated mRNA species given the primers used in each assay is given in green. Products above this length originate from 
polyadenylated mRNA molecules. Marked differences in poly(A) tail length distribution are detected for Smed-agat-1 , showing an increased frequency 
of long poly(A) tails after Smed-not1 knock down. Slight differences are observed for Smed-nb.21.11e and Smed-eif-3. 
doi:10.1371/journal.pgen.1004003.g005 



that neoblast maintenance is not affected by Smed-notl knock down 
at these time points. 

Since WMISH does not provide a quantitative measure of 
mRNA levels we quantified these differences by qRT-PCR 
experiments on RNA from Smed-notl (RNAi) animals. Smedtud-1, 
Smed-vasa-1 and Smed-pcna all progressively increased to levels 
approximately 50% and 100% higher than those of control (RNAi) 
after 10 and 15 days, respectively (Figure 6D). Smedwi-1 
transcript levels were also significantly increased in whole 
animals after 15 days (Figure 6D). Collectively, these results 
demonstrate that Smed-notl RNAi knock down leads to an 



increase in mRNA levels in several genes expressed in neoblasts 
and their progeny. 

Given the known conserved function of the CCR4-NOT 
complex in regulating mRNA levels through targeted dead- 
enylation we performed PAT assays on the set of neoblast 
markers and on the samples above (Figure 5D). In all cases Smed- 
notl knock down resulted in increased average poly(A) tail length 
(Figure 6E), demonstrating that increased transcript levels 
correlate with a failure in deadenylation. These data confirm 
that knock down of the CCR4-NOT complex subunit Smed-notl 
leads to increased transcript levels of genes known to be key to 
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Figure 6. Smed-notl(RNAi) animals have increased levels of neoblast transcripts with increased frequency of long poly(A) tails. (A-l) 
WMISH of the neoblast marker Smedtud-1 (A), Smed-vasa-1 (B) and Smed-pcna (C) in control(RNAi) animals (upper panels) and Smed-notl(RNAi) 
animals 1 0 (middle panels) and 1 5 (bottom panels) days after RNAi, showing normal expression of these mRNAs after Smed-notl knock down, though 
qualitative differences in the level of expression are suggested. Anterior is to the left. Scale bars: 500 ^m. (D) Quantification of the level of expression 
by qRT-PCR of the neoblast markers Smedwi-1, Smedtud-1, Smed-vasa-1 and Smed-pcna in control(RNAi) (c) and Smed-notl(RNAi) animals (n) 10 and 15 
days after RNAi, normalized expression and relative to respective control(RNAi) samples. Error bars represent standard deviation, asterisks represent 
statistical significance. Smedtud-1, Smed-vasa-1 and Smed-pcna transcripts accumulate progressively after 10 and 15 days of RNAi, while Smedwi-1 
only accumulates significantly after 15 days. (E) PAT assays reflecting the distribution of mRNA poly(A) tail lengths for Smedwi-1, Smedtud-1, Smed- 
vasa-1 and Smed-pcna in control(RNAi) (c) and Smed-notl (RNAi) animals (n) 10 and 15 days after RNAi. Size markers used are represented in blue, the 
theoretical length of the amplicon corresponding to the deadenylated mRNA species given the primers used in each assay is given in green. Marked 
differences in poly(A) tail length distribution are detected for all four mRNAs, showing an increased frequency of long poly(A) tails after Smed-notl 
knock down. 

doi:10.1371/journal.pgen.1004003.g006 



neoblast function, by blocking their deadenylation and subse- 
quent degradation. 

Upregulation and increased polyadenylation of Smedtud- 
7, Smed-vasa-1 and Smed-pcna occurs specifically in stem 
cells 

Since many neoblast mRNAs are also expressed in differenti- 
ated cells (e. g. Smedtud-1 and Smed-vasa-1 are also prominently 
expressed in the CNS, Figure IB— C), the increased mRNA levels 
we detected could arise from a response in differentiated cells 
alone, stem cells alone or from both differentiated and stem cells. 
To distinguish these possibilities, we planned to use an irradiation 



approach to remove all neoblasts and then measure transcript 
levels in Smed-notl (RNAi) and control(RNAi) worms by qRT-PCR. 
We reasoned that if transcript accumulation was indeed limited to 
stem cells then mRNA levels of these transcripts after irradiation 
would be equal in both irradiated control(RNAi) and irradiated 
Smed-notl (RNAi) samples. 

Smedtud-1 and Smed-vasa-1 are expressed in the CNS to levels 
that amount respectively to roughly 70% and 40% of their total 
expression, according to our previous neoblast profiling by a 
combinatorial RNA-seq, RNA interference and irradiation 
approach [50]. We confirmed this by qRT-PCRs (Figure S5A) 
and observed that Smed-pcna expression amounts to roughly only 
10% of its normal expression 24 hours after irradiation (Figure 
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S6A), indicating that 24 hours of irradiation suffice to eliminate 
around 90% of neoblasts. Conversely, most neoblast progeny 
survive beyond 1 day post-irradiation, as measured by our qRT- 
PCR experiments with the markers Smed-nb.21.1 le and Smed-agatl 
(Figure S6B) and consistent with previously published data 
[50,66], making then 24 hours after irradiation the ideal time 
point to perform our experiment. 

Therefore, in order to find out if the increased levels of neoblast 
mRNAs originate in neoblasts or instead in the CNS or elsewhere, 
we then used this irradiation approach in Smednot- 1 (RNAi) animals 
and compared them to control(KNAi) animals. We irradiated both 
Smednot- 1 (RNAi) and control(RNAi) animals at 9 and 14 days after 
dsRNA administration, 24 hours before the data collection time 
points of 10 and 15 days respectively. We then examined the 
expression pattern of Smedtud-1 by WMISH. This confirmed our 
qRT-PCR experiments, and was consistent with our previous 
WMISH results (Figure 7A). Non irradiated Smed-notl (RNAi) 
animals showed qualitatively more intense expression of Smedtud-1, 
However, all animals that were irradiated 24 hours prior to 
fixation showed an identical expression pattern of Smedtud-1, with 
similar levels of signal detected only in the CNS to their 
control(RNAi) irradiated counterparts (Figure 7A). This experiment 
shows that Smedtud-1 is not ectopically expressed in other tissues or 
organs after Smed-notl knock down, since this ectopic expression 
should be visible either in non irradiated or irradiated samples, 
and suggests instead that the increased levels of Smedtud-1 come 
from an accumulation of this transcript in neoblasts. 

We performed qRT-PCR and confirmed that Smedtud-1, Smed- 
vasa-1 and Smed-pcna mRNA levels increased progressively in Smed- 
notl (RNAi) animals (Figure 7B, left) but both control(RNAi) and 
Smed-notl (RNAi) animals irradiated 24 hours previously contained 
similar levels of all three transcripts (Figure 7B, right). Similar to 
results for wild type worms (Figure S6A), the levels of the Smedtud- 
1, Smed-vasa-1 and Smed-pcna transcripts were respectively around 
70%, 40% and 10% of the expression in non-irradiated 
control(RNAi) animals. Therefore we conclude that the overexpres- 
sion observed in Smed-notl (RNAi) animals for these three 
transcripts disappears 24 hours after irradiation, and is therefore 
located in irradiation-sensitive neoblasts, rather than elsewhere in 
the body. 

To implicate the CCR4-NOT complex-mediated targeted 
mRNA degradation directly in this effect we performed PAT 
assays on non irradiated and irradiated samples and observed that 
the distribution of poly(A) tail lengths of all these neoblast mRNAs 
was increased in Smed-notl (RNAi) animals, but this effect was 
removed by irradiation (Figure 7C and D). In the case of neoblast 
specific Smedwi-1 and Smed-pcna no signal was detected after 
irradiation (Figure 7C). In the case of Smed-vasa-1 and Smedtud-1 
signal was detectable after irradiation, reflecting expression in the 
CNS, but these transcripts did not show any increase in poly(A) tail 
length distribution after Smed-notl knock down (Figure 7D). These 
data confirm that effects on mRNA levels for these important 
neoblast genes are confined to stem cells. As an additional control 
we also measured Smed-nb.21.1 le and Smed-agat-1 poly(A) lengths 
(Figure 7E). 24 hours after irradiation cells expressing these 
transcripts are still present and the increase in poly(A) tail length 
caused by Smed-notl knock down is still evident. This confirms that 
irradiation itself does not cause the absence of detected differences 
in poly(A) tail length per se, but by removing the cycling neoblasts. 
The poly(A) tail length distribution of control mRNAs were not 
affected by irradiation (Figure 7F). 

As a further proof that the effects we observe are confined to 
stem cells we looked at gene expression, poly(A) tail length and the 
efficacy of Smed-notl knock down across FACS cellular compart- 



ments [5 1,72]. Progressive increases in neoblast gene mRNA levels 
were confined to XI and X2 populations of sorted cells from Smed- 
notl (RNAi) animals (Figure 8A, XI and X2). Both of these sorted 
fractions contain neoblasts to different extents. Levels of these 
mRNAs were not increased in irradiation insensitive (Xins) sorted 
fractions from Smed-notl (RNAi) animals compared to controls 
(Figure 8A, Xins). This fraction contains primarily differentiated 
cells including CNS cells which also express Smedtud-1 and Smed- 
vasa-1. 

One possibility for the specificity we observe is that Smed-notl 
knock down has a high efficacy in stem cells but not in post-mitotic 
cells, as it has been suggested recendy for Smed-bruno-like knock 
down [38,73]. We performed qRT-PCR measurement of Smed- 
notl transcript levels and found that Smed-notl mRNA is 
consistently expressed in all cellular fractions (consistent with 
Figure S2 and Figure 1A) and that Smed-notl knock down 
significandy depletes Smed-notl in all sorted compartments 
compared to control(RNAi) animals (Figure 8B). Consequendy the 
observed neoblast specificity is unlikely due to an absence of knock 
down in differentiated cells. 

To finally link the mechanism of CCR4-NOT complex- 
mediated deadenylation to increased mRNA levels we also 
checked poly(A) tail lengths in FACS sorted cells. For Smedtud-1, 
Smed-vasa-1, Smed-pcna and Smedwi-1 we observed a progressive 
increase in long poly(A) tails in the XI compartment of stem cells 
(Figure 8C and D). This trend was also observed for these 
transcripts in X2 cells (Figure 8C and D). No poly(A) tail signal 
was detectable in Xins fractions for Smed-pcna and Smedwi-1, 
consistent with their low abundance in this fraction, and poly(A) 
tail length was not affected for Smedtud-1 and Smed-vasa-1 in Xins 
cells (Figure 8C and D). These data are in agreement with our 
irradiation based experiments (Figure 7). The poly(A) tail lengths 
of control mRNAs Smed-ef-2 and Smed-mhc were only mildly or not 
affected by Smed-notl knock down (Figure 8E). These mild 
differences are likely due to CCR4-NOT mediated deadenylation, 
but again seem to be restricted to XI and X2 cells, whereas Xins 
cells remain unaffected. To which extent all, most or only a subset 
of neoblast transcripts are affected after Smed-notl knock down 
remains an open question. 

Taken together, these results clearly demonstrate that Smed-notl 
knock down induces a prominent increase of key transcripts 
expressed in neoblasts, that this accumulation occurs specifically in 
neoblasts, and that it is associated with an increased frequency of 
long poly(A) tails of these transcripts specifically in neoblasts. 
Furthermore, these results strongly suggest that the neoblast- 
specific increase of mRNA levels of genes such as Smedtud-1, Smed- 
vasa-1 and Smed-pcna may be responsible for the impaired 
differentiation capacities of neoblasts observed in Smed-notl (RNAi) 
animals. It is likely that other genes expressed in neoblasts are 
similarly upregulated in Smed-notl (RNAi) animals and contribute to 
differentiation defects. Our results suggest a mechanism by which 
the differentiation capacities of neoblasts are dependent on CCR4- 
NOT mediated degradation of specific neoblast mRNAs. 

Discussion 

Planarians are an emerging in vivo model for stem cell biology 
because of their unique stem cell population. In this study we used 
the planarian S. meditenanea as a model system to establish a 
function for the CCR4-NOT complex in stem cell regulation. We 
found that Smed-notl knock down abrogated regeneration and 
impaired homeostatic cell turnover. Interestingly, Smed-notl knock 
down primarily affects the stem cell compartment of S. meditenanea 
rather than inducing more widespread effects, even though the 
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Figure 7. Increased levels of neoblast transcripts and their increased frequency of long poly(A) tails are irradiation sensitive. (A) 

WMISH of the neoblast marker Smedtud-1 in control(RNAi) animals (left panels) and Smed-notl(RNAi) animals (right panels) non irradiated (top panels) 
and irradiated 24 hours (bottom panels) prior to data collection time points 1 0 and 1 5 days after RNAi. A consistent qualitative difference is observed 
in non irradiated animals, however, no qualitative differences are observed in irradiated animals, showing that Smed-notl (RNAi) animals do not 
overexpress ectopically Smedtud-1. Anterior is to the left. Scale bars: 500 jj.m. (B) Quantification of the level of expression by qRT-PCR of the neoblast 
markers Smedtud-1, Smed-vasa-1 and Smed-pcna in control(RNAi) (c) and Smed-not1 (RNAi) (n) animals non-irradiated and irradiated 24 hours prior to 
data collection time points 10 and 15 days after RNAi, normalized expression and relative to respective control(RNAi) samples. Error bars represent 
standard deviation, asterisks represent statistical significance. Smedtud-1, Smed-vasa-1 and Smed-pcna transcripts accumulate progressively after 10 
and 1 5 days of RNAi, but this accumulation is eliminated 24 hours post irradiation, with Smed-notl (RNAi) irradiated animals showing levels similar to 
control(RNAi) irradiated animals. (C-F) PAT assays reflecting the distribution of mRNA poly(A) tail lengths for the neoblast specific mRNAs Smedwi-1 
and Smed-pcna (C), the neoblast and CNS expressed mRNAs Smed-vasa-1 and Smedtud-1 (D), the progeny specific mRNA Smed-nb.21.1 1e and Smed- 
agat-1 (E) and the housekeeping and tissue specific mRNAs Smed-ef-2, Smed-eif-3, Smed-mhc and a spiked-in control in control(RNAi) (c) and Smed- 
notl (RNAi) (n) animals non irradiated and irradiated 24 hours prior to data collection time points 15 days after RNAi. Size markers used are 
represented in blue, the theoretical length of the amplicon corresponding to the deadenylated mRNA species given the primers used in each assay is 
given in green. Neoblast specific markers are not detected after irradiation (C). The marked differences in poly(A) tail length distribution detected for 
the neoblast and CNS mRNAs Smed-vasa-1 and Smedtud-1 are eliminated by irradiation, showing that the fractions of these mRNA populations 
localised in the CNS show no differences after Smed-notl knock down (D). The differences in poly(A) tail length distribution detected for the progeny 
specific mRNAs Smed-nb.21.1 1e and Smed-agat-1 are not affected by irradiation, as these cells are not eliminated after 24 hours of irradiation (E). No 
differences are detected for Smed-ef-2, Smed-eif-3, Smed-mhc and a spike-in control RNA. 
doi:1 0.1 371 /journal.pgen.1 004003.g007 



CCR4-NOT complex is the major deadenylating complex in 
eukaryotes and regulates at least 85 % of mRNAs in yeast [74] . 

While we observed a stark and specific effect of Smed-notl knock 
down on deadenylation, it is still possible that other functions of 
the CCR4-NOT complex might also be impaired. The CCR4- 
NOT complex is involved in several steps of RNA metabolism 
[12,13,16] and further work is therefore needed to elucidate which 
ones are also at work in stem cells. 

We observed effective gene knock down of Smed-notl even in 
differentiated cell fractions, but specific effects on neoblast 



transcripts were limited to stem cells. This fact suggests that 
targeted deadenylation by either RNA-binding proteins or 
miRNAs is providing specificity and is therefore central to stem 
cell differentiation and self-renewal properties. Consistently, 
several studies have implicated the CCR4-NOT complex in 
mRNA-specific deadenylation via targeted recruitment of the 
CCR4-NOT complex by RNA-binding proteins [26,27,30,31], 
which are in turn known to be highly enriched and functionally 
important in neoblasts [43-46,48-50,52,75]. It is possible to 
hypothesize that disruption of the CCR4-NOT complex via knock 
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Figure 8. Increased level of transcripts and increased frequen- 
cy of long poly(A) tails are restricted to neoblast-containing 
cell populations. (A-B) Quantification of the level of expression by 
qRT-PCR of the neoblast markers Smedwi-1, Smedtud-1, Smed-vasa-1 
and Smed-pcna (A) and Smed-notl (B) in FACS sorted populations X1, X2 
and Xins in control(RNAi) (c) and Smed-notl(RNAi) (n) animals 10 and 15 
days after RNAi, normalized expression and relative to respective XI 
control(RNAi) samples. (A) Smedwi-1, Smedtud-1, Smed-vasa-1 and Smed- 
pcna transcripts accumulate progressively after 10 and 15 days of RNAi 
in X1 and X2 cells, the two fractions that contain neoblasts to different 



extents, but this accumulation is not observed in Xins cells, which 
contain differentiated cells exclusively, including CNS cells. (B) Smed- 
notl is significantly depleted across all three cell fractions, showing that 
the absence of accumulation and increased frequency of long poly(A) 
tails of neoblast mRNAs that are expressed also in CNS is not due to 
absence of effective gene knock down in differentiated cells. Error bars 
represent standard deviation and asterisks represent statistical signif- 
icance in A-B (C-E) PAT assays reflecting the distribution of mRNA 
poly(A) tail lengths for the neoblast and CNS expressed mRNAs 
Smedtud-1 and Smed-vasa-1 (C), the neoblast specific mRNAs Smed- 
pcna and Smedwi-1 (D) and the housekeeping and tissue specific 
mRNAs Smed-mhc and Smed-ef-2 (E) in FACS sorted populations XI, X2 
and Xins from control(RNAi) (c) and Smed-notl(RNAi) (n) animals 10 and 
15 days after RNAi. Size markers used are represented in blue, the 
theoretical length of the amplicon corresponding to the deadenylated 
mRNA species given the primers used in each assay is given in green. 
(C). The marked differences in poly(A) tail length distribution detected 
for the neoblast and CNS mRNAs Smed-vasa-1 and Smedtud-1 are only 
detected in XI and X2 but not in Xins FACS sorted populations, 
showing that the fractions of these mRNA populations localised in the 
CNS show no differences after Smed-notl knock down. (D) The marked 
differences in poly(A) tail length distribution detected for the neoblast 
specific mRNAs Smed-pcna and Smedwi-1 are only detected in X1 and 
X2 but the mRNAs are not detected in Xins FACS sorted populations. (E) 
No differences in poly(A) tail length distribution are detected for the 
tissue specific mRNA Smed-mhc, and only slight differences are 
detected in XI and X2 but not in Xins fractions for the housekeeping 
mRNA Smed-ef-2. 

doi:10.1371/journal.pgen.1004003.g008 

down of its scaffolding subunit might impair the protein-protein 
interactions needed to tether deadenylating activity to tran- 
scripts targeted by RNA-binding proteins or miRNAs for 
degradation, while general non-targeted deadenylation would 
remain relatively unaffected, therefore not causing a broader 
metabolic failure at earlier time points. In fact, the miRNA 
silencing complex protein GW182 specifically interacts with the 
Notl subunit of the CCR4-NOT complex through specific and 
conserved domains [36,37]. 

Despite the specificity of the effects seen in neoblast transcripts, 
mRNAs expressed elsewhere were also found to be affected. Both 
Smed-agat-1 and, to a lesser extent, Smed-nb.21.11e were also 
affected. These results show that CCR4-NOT deadenylating 
activity is present in cell types other than neoblasts and that 
specificity is not due to restriction of activity to stem cells. 
Furthermore, the effects on well described neoblast progeny 
markers suggest that Smed-notl knock down likely influences several 
steps of cellular differentiation that may all contribute to the 
observed effects on homeostasis and regeneration. 

The Smed-notl phenotype is progressive with respect to both the 
decreasing capacity of the animals to produce blastema cells and 
by the accumulation of mRNAs in stem cells. The phenotype 
results from a drop in neoblast progeny numbers, followed by stem 
cell loss. Similar neoblast and progeny dynamics have been shown 
by us in Smedwi-2 (RNAi) organisms here and by another group in 
Smed-CHD4(RNAi) organisms [64]. Both Smedwi-2 and Smed- 
CHD4 knock downs initially affect neoblast differentiation rather 
than their self-renewal and proliferative capabilities [44,64] which 
are only affected at later time points, similar to Smed-notl knock 
down. The ultimate cause for stem cell loss after the impairment of 
neoblast differentiation is unknown, and likely to be a broad failure 
in homeostasis as organs and tissues fail. However, for more than 
15 days after Smed-notl dsRNA administration proliferating 
neoblasts are detectable in large numbers while regeneration is 
abrogated, suggesting neoblast loss is not a primary cause for the 
regeneration defect, instead implicating impaired neoblast differ- 
entiation capacities. 
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Smed-notl knock down induces an increased frequency of long 
poly(A) tails of Smedtud-1, Smed-vasa-1, Smed-pcm and Smedwi-1 
mRNAs as these same mRNAs accumulate. This regulation occurs 
specifically in neoblasts, rather than in the CNS, where two of 
these transcripts are also expressed. However, it is likely that Smed- 
notl knock down affects many more transcripts that contribute to 
failure in stem cell differentiation. Due to the lack of specific 
antibodies it is difficult to evaluate if the differences observed in 
transcript abundance and polyadenylation state affect the abun- 
dance of the proteins that these transcripts encode. However, 
taking into account that all four transcripts expressed in neoblasts 
analyzed accumulate in these stem cells, it is reasonable to 
conclude that Smed-notl knock down induces transcriptome-wide 
changes in stem cell expression patterns, and that these changes 
will likely affect protein levels. 

Our results offer a new mechanistic insight into post-transcrip- 
tional regulation in neoblasts and its targets. After depletion of a 
post-transcriptional regulator many transcripts accumulate with- 
out being degraded, and this likely prevents neoblast differentia- 
tion, which needs the effective removal of these transcripts. The 
deadenylating activity of the CCR4-NOT complex is clearly 
central to this process. For genes like Smedtud-1 and Smed-vasa-1 , 
which are expressed in neoblasts and the CNS, we observe 
accumulation and increased frequency of long poly(A) tails of the 
transcript only in neoblasts. This suggests that the deadenylation of 
these mRNAs is regulated specifically during the onset of 
differentiation and requires the targeted recruitment of the 
CCR4-NOT complex by RNA-binding proteins, as has been 
described in other organisms [26,27,30,31,76]. Interestingly, 
several RNA binding proteins and post-transcriptional regulators 
have already been described as crucial for neoblast function 
[43,44,46,48], and some of them have been already been 
functionally linked with the CCR4-NOT complex in other model 
organisms. Future research will help elucidate the mechanisms by 
which these proteins orchestrate planarian stem cell processes. 

The CCR4-NOT complex has been shown to mediate the 
repressive function of both miRNAs and piRNAs [32-35]. Small 
RNAs are believed to be very important regulators of mammalian 
stem cells [77] and neoblasts [78,79]. The Smed-notl RNAi 
phenotype is very similar to those of Smedwi-2 and Smedwi-3, two 
Piwi proteins involved in piRNA regulation [44,80] . Furthermore, 
several studies have highlighted the presence of miRNAs highly 
enriched in stem cells [78,79,81]. Future research will help in 
understanding if these phenotypic similarities reflect a functional 
link between Piwi proteins, piRNAs, miRNAs and the CCR4-NOT 
complex in planarian stem cells. Our results highlight the 
importance of the CCR4-NOT complex in the regulation of stem 
cells, the fact that post-transcriptional regulation of gene expression 
is a key element in the regulation of pluripotency, and that 
planarians will provide an excellent platform for these studies. 

Materials and Methods 

Organisms 

Planarians of the asexual strain of S. mediterranea were kept and 
used as previously described [82] . 

Sequences 

The putative members of the S. mediterranea CCR4-NOT 
complex and other transcripts were identified in published S. 
mediterranea transcriptomic and genomic sequences [53,58-60] and 
the longest transcripts for Smed-not-1 confirmed by PCR and 
sequencing. The putative members of the S. mediterranea CCR4- 
NOT complex were identified by TBLASTN searches in the 



current assembly of the S. mediterranea genome and in the available 
transcriptomic data. In order to determine the number of loci for 
each of the components, the different transcripts identified were 
mapped to the S. mediterranea genome. The genomic sequence 
encoding Smed-notl was found split in two contigs (v3 1.00 1778 and 
v3 1.002774), the existence of one single transcript for these two 
genomic contigs was confirmed by PCR using the primers 5'- 
CATCGCAACAATGGAGAGAA-3' and 5'-ATTTGAGCTG- 
TATGGGCGAT-3', each mapping to one of the two contigs. 
These PCR experiments revealed as well the existence of a 3 kb 
region not present in the S. mediterranea genomic data. 

The full sequence of Smed-notl was obtained by de novo 
assembling the transcript from the raw transcriptomic data, using 
the known transcriptomic and genomic data as a guide. The Smed- 
notl sequence has been deposited in Genbank (accession 
KF781122). 

The sequences of Smedtud-1, Smed-vasa-1 (accession JQ425 140) 
and Smed-pcna (accession EU856391) were found in our reference 
transcriptomic data, encoded by the transcripts AAA.454ES- 
TABI.16133, AAA.454ESTABI. 18605 and AAA.454ES- 
TABI.22122 respectively. Smed-ef-2 is encoded by the transcript 
AAA.454ESTABI. 17328. The Smedtud-1 sequence has been 
deposited in Genbank (accession KF781126). 

RNAi 

RNAi experiments were performed as previously described 
[82]. Control(RNAi) worms were injected with dsRNA encoding for 
GFP, a gene not present in the S. mediterranea genome. dsRNA 
encoding for Smed-notl was prepared by in vitro transcription of a 
region of the Smed-notl gene. Briefly, an amplicon was generated 
from S. mediterranea reverse transcribed RNA with the primers 5'- 
GGCCGCGGTGTCCAAGAAAAAGC AAGTCAG-3 ' and 5'- 
GCCCCGGCCAGCTGGCGTCAGTTTAGTGAA-3', contain- 
ing a 5' adaptor sequence for the subsequent addition of T7 
promoter. The product of this PCR was gel purified and subjected 
to another step of amplification with the primers 5'-GAGAATTC- 
TAATACGACTCACTATAGGGCCGCGG-3 ' and 5'-AGG- 
GATCCTAATACGACTCACTATAGGCCCCGGC-3' with 
the purpose of attaching T7 promoter sequences to both ends of 
the amplicon. The product of this PCR was further purified with 
Purelink DNA purification columns (Invitrogen) and used as a 
template for in vitro transcription using T7 RNA polymerase 
(Roche). The product of the in vitro transcription was treated with 
Turbo DNAse (Ambion), phenol extracted, precipitated in ethanol 
in presence of sodium acetate, glycogen and EDTA and re- 
suspended in water. dsRNA encoding GFP for use as a negative 
control was similarly prepared from a vector encoding the GFP 
gene. The final concentration of the injected solution was 1 p.g/ 
(Xl. Animals were injected with a Nanoject II (Drummond) for 
three consecutive days and monitored or used for experiments in 
the subsequent days. Day 1 after RNAi is considered to be in all 
experiments the first day after the third dsRNA injection. Alive 
animals were imaged in a Zeiss Discovery V8 with a Zeiss 
AxioCam MRC camera. 

Irradiation 

Irradiation was performed as previously described [79]. Animals 
were placed in a sealed y-ray source and administered an 
irradiation dose of 100 Gy. 

In situ hybridization, immunohistochemistry and imaging 

WMISH, FWMISH and IHC were performed and imaged as 
previously described [82] . When qualitative differences are shown, 
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animals were processed and monitored in parallel. Riboprobes 
were generated by in vitro transcription of PCR products 
generated as described above, with only one T7 promoter linked 
to the 3' end of the amplicon, and in the presence of digoxigenin- 
labelled UTP (Roche). The products of in vitro transcription 
reactions were then treated with Turbo DNAse (Ambion). 
Riboprobes were then precipitated in ethanol in the presence of 
LiCl and glycogen and resuspended in 50% formamide in TE 
buffer, 0.01% Tween. 

The following primers were used: 

Srmd-noth 5'- GGCCGCGGTGTCCAAGAAAAAGCAAGT- 
CAG-3' 

and 5'- GCCCCGGCCAGCTGGCGTCAGTTTAGTGAA-3' 
Smedtud-1: 5'-GGCCGCGGCTAATGCCAGTTGACTGTC- 
C-3' 

and 5'-GCCCCGGCCCGAAAAAGTTCCGCATCACTT-3' 
Smedvas-L 5 ' -GGC C GCGGAGCTGTTGGAGTTGTTGG- 
CTCAG-3' 

and 5'-GCCCCGGCCCTAATCTTCGAGCCATTCAG-3' 
Smed-pcna: 5 ' -GGC C GCGGATGGACTTGGATGGAGAT- 
CACT 

Smedwi-1: 5 ' -GGC C GCGGAAGTGGTGGTATTCGAGAA- 
GGA-3' 

and 5'-GCCCCGGC C ACGAATC GTAATC GGTTGTTCT-3 ' 
Smed-agat-1: 5'-GGCCGCGGGAAATGATTGAGTCCACC- 
ATGA-3' 

and 5 '-GCCCCGGCCTGCAATATCTGGATAAGGAGCA-3 ' 
Smed-nb.21.Ue: 5'-GGCCGCGGGTGATTGCGTTCGCG- 
T AT ATT- 3' 

and 5 '-GCCCCGGCCATTTATCCAGCGCGTCATATTC-3 ' 
Briefly, animals were killed in a 2% HQ solution, fixed in 
Carnoy's solution, bleached in a 8% H 2 0 2 /methanol, rehydrated, 
permeabilized with Proteinase K (Sigma), treated with 0.25% and 
0.5% acetic anhydride in 0.1M triethanolamine pH 7.6, prehy- 
bridized and hybridized with digoxigenin labelled riboprobes 
(0.2 ng/ul, O/N at 56°C) They were then washed in buffers of 
increasing stringency, immunolabelled with anti-digoxigenin- 
alkaline-phosphatase antibody (Roche) and developed in the 
presence of NBT and BCIP (Roche). For FWMISH, an anti- 
digoxigenin-peroxidase antibody (Roche) was used and the signal 
was developed with the Tyramide Signal Amplification kit (Perkin 
Elmer). WMISH specimens were imaged on a Zeiss Discovery V8 
equipped with a Zeiss AxioCam MRC camera. FWMISH 
specimens were imaged in a Leica SP3 confocal. 

Quantification of Smed-agat- 7 -positive cells 

Animals (N = 8 per time point and treatment) were stained by 
FWMISH of Smed-agat-1 and monitored by confocal microscopy. 
Cell counts were performed in z-projections of both the dorsal and 
ventral sides of the animals. 130 squares corresponding to 
250 um 2 in both the dorsal and ventral parts of the animals were 
selected for counting along the length of the animals. Counts were 
performed using ImageJ software. Error bars represent standard 
error of the mean. Statistical significance was analyzed by 
Student's T test by comparing values from each sample to its 
respective control sample. 

Immunohistochemistry 

Whole mount immunohistochemistry of phospho-histone-3 was 
performed as previously described [83]. Essentially, animals were 
fixed as above, blocked in a 1% BSA/PBS 0.3% triton X-100 
solution, incubated overnight with anti-phospho-histone-3 (Milli- 
pore, 1/500 dilution) and anti-phospho-tyrosine (Cell Signaling, 
1/200 dilution), washed and incubated in Alexa Fluor 488 and 



568 secondary antibodies (Molecular Probes, 1/400 dilution). 
Animals were then washed, mounted in 70% glycerol/PBS and 
imaged with a Leica SP3 confocal and a Leica MZ16F 
fluorescence stereomicroscope and a Leica DFC 300Fx camera. 
Countings of phospho-histone-3 were performed with ImageJ 
software and normalized to the total area of the sample. 

qRT-PCR 

qRT-PCR experiments were performed as previously described 
[48] with modifications. Essentially, total RNA from samples of 5 
animals was extracted with Trizol reagent (Invitrogen) according 
to the manufacturer's instructions, and cDNAs were synthesized 
with SuperScriptHI Reverse Transcriptase (Invitrogen). qRT-PCR 
experiments were then performed using the Absolute qPCR SYBR 
Green Master Mix (Thermo Scientific). Experiments were 
performed on three biological replicates per time point and 
treatment. Each biological replicate was technically replicated 
three times in each reaction, each reaction was repeated three 
times. The gene Smed-ef-2 was used for normalization. 

The following primers were used: 

Smed-notl: 5 ' -GAC AGC GATTATGAACTGCC -3 ' and 5'- 
CTGCTATGTTACTGGTGTTGAG-3' 

Smedwi-1: 5'-AGTTCCTGTTCCAACGCATTATG-3' and 
5 '-CTGGAGGAGTAACACCACGATGA-3 ' 

Smed-nb.21.11e: 5'-GTCTCCCGCCAAATCAAGTA-3' and 
5 '-TTTCATGCAATCTGCTTTCG-3 ' 

Smed-agat- 1:5' -TC CATC C AG AAC C G ATTG AT- 3 ' and 5'- 
CTC C C AAGTC ATGGTGG AC T-3 ' 

Smedtud-1: 5 ' -TGATGAAGGAACTTC GGGTGAT-3 ' and 5'- 
TCTGAGCAACCGATTGAAACC-3' 

Smedvas-1: 5 ' -TGAAATGAAC AAATC C C GAC-3 ' and 5'-GA- 
GAGC C AAAC T AATTC C AG- 3 ' 

SWrf-/>OTa:5'-GGCGCTTGGTAGTAATGATTCCCTA-3' 
and 5'-TACCTAAGTGATCTCCATCCAAGTCC-3' 

Sm(;rf-ef-2:5'-CAGCCAGTAGCTTTAAGCGATGA-3' and 5'- 
ACTCTCAACGCTGCTGTC ACTTC-3 ' 

Statistical significance was measured by Student's T test by 
comparing values from each sample to their respective control 
sample. 

FACS 

Fluorescence activated cell sorting (FACS) of planarian samples 
was performed as previously described [51,72]. Analysis was 
performed with Flowjo. 

PAT assays 

PAT assays were performed as previously described [69,71] 
with minor modifications. A total of 400-1000 ng of total RNA 
extracted from five animals were used per each time point, 
replicate and treatment, except for FACS samples, in which 40 
animals were used per time point and treatment and PAT 
reactions were performed with 100 ng of total RNA. Three 
biological replicates were analyzed per time point and treatment, 
and technically replicated at least twice, except for FACS samples, 
which were only technically replicated. C. elegans total RNA was 
spiked- in as a control and assayed with a primer for cpg-2. 

RNA samples were incubated with 0.3 ug of 5'-phosphorylated 
oligo d(T) in a total volume of 8 ul, and heat denatured for 5 min 
at 65°C. Then, the following mixture was added: 4 ul of Super 
Script II First Strand Buffer, 0.5 ul of 0.1 M DTT, 2.25 ul of 
10 mM ATP, 0.125 ul of RNAsin (Promega), 1.25 ul of 1 mM 
(each) dNTPs (Promega), and 1 ul of 2000 units/ul T4 DNA 
Ligase (New England Biolabs). The volume was then brought to a 
total of 20 ul with water, and the samples were incubated at 42°C 
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for 30 min. in order to allow the saturation of polyA tails with 
oligo d(T). Then, 0.5 (0.1 of an oligod(T)-anchor 100 uM was added 

the mixture was incubated for 2 hours at 12°C to allow ligation of 
oligo d(T) molecules. Then, samples were pre-warmed at 42°C 
and 1 ul of Superscript II (Invitrogen) was added. Finally the 
samples were incubated for 1 hour at 42°C and the Superscript 
enzyme was heat inactivated at 70°C for 30 min. 

These cDNAs were then used in PCR reactions, with fresh 
aliquots of the anchor primer and gene specific forward primers 
designed close to the 3 ' region of the mRNA. The sequences of the 
mRNAs tested were obtained in our reference transcriptome (4). 
The following primers were used: 

Smedtud-h 5'-TGATGAAATAATGCTACCCGCGCAAT-3' 
Smed-vasa-1: 5'-AGCCGACTTCTGAATGGCTCGAAGA-3' 
Smed-pcna: 5 '-CAAAGGCTGC ACCTCTTTCTTCTC A-3 ' 
Smedwi-1 : 5 ' -C GTTGGC AAGATTC ATC GTGGTGTT-3 ' 
Smed-agat-1: 5'-TCGGATGTTAGAAGGCGAGGAGACC-3' 
Smed-nb.21.Ue: 5'- GACGGCCAGTGTGACGCAGAAT-3' 
Smed-ef-2: 5'-AACCCACTGGATCCCACAACGAAAC-3' 
Smed-eif-3: 5 ' -GTTGCCC CATC GATTGGATACTTC G-3 ' 
Smed-mhc: 5'- CGAGGAGCAAGTTCTGGACCTGGAA-3' 
PCR reactions were carried away with DreamTaq (Fermentas) 
and amplified for 28-32 cycles of 94°C for 20 s., 65°C for 20 s. 
and 72 for 30 s. The products of the PCR reactions were analyzed 
on 1.5% agarose gels. 

Supporting Information 

Figure SI CCR4-NOT complex deadenylases do not induce a 
strong phenotype in planarians. (A) WMISH of Smed-not6, Smed- 
not7A and Smed-not7B in non irradiated and 3 days post irradiation 
animals. (B) control(KNAi) , Smed-not6(RNAi), Smed-not7A(KNAi) and 
Smed-not7B(RNAi) animals cut 5 days after RNAi and monitored 6, 
9 and 19 days after RNAi. Only a weak phenotype of delayed 
regeneration is observed for Smed-not7A. (C) WMISH of Smed-H2B, 
Smed-nb.21.11e and Smed-agat-1 in amtrol(KMAi) , Smed-notG(RNAi), 
Smed-not7A(RMAi) and Smed-not7B(KNAi) animals 10 days after 
RNAi. No alteration of neoblast or progeny markers is observed. 
(TIF) 

Figure S2 Smed-notl is highly expressed across all planarian 
FACS sorted fractions. Expression levels, XI vs Xins enrichment 
and corresponding gene IDs of the transcripts encoding for 
Smedtud-1, Smed-vasa-1, Smed-pcna, Smed-eJ2 and Smed-notl. Data 
taken from Onal et al. 2012. When transcripts are split into 
different transcriptomic sequences, all sequences are shown 
independendy. The neoblast expressed transcripts Smedtud-1, 
Smed-vasa-1 and Smed-pcna are most highly expressed in XI 
fractions. The enrichment vs. Xins fractions is most high in Smed- 
pcna and lower in Smed-vasa-1 and Smedtud-1, consistent with their 
expression in CNS. The enrichment vs. Xins fractions of Smed-notl 
is lower than all three neoblast expressed transcripts and more 
similar to the housekeeping gene Smed-ef-2. XI /Xins: log2(RPKM 
Xl)-log2(RPKM Xins). 
(TIF) 

Figure S3 Smed-notl is required for planarian regeneration and 
homeostatic cell turnover. (A-F) Control(KNAi) (A) and Smed- 
notl (RNAi) animals cut 1 (B), 3 (C), 5 (D), 10 (E) and 15 (F) days 
after RNAi, and monitored every 2 days after transection. All 
panels are anterior wounds. Time of regeneration is indicated on 
top, total days after RNAi are indicated in each panel. Five 
animals were used per time point. 5 control(RNAi) animals were 
used for each of the time points, only one is shown (1 day) since no 



differences were detected among them. Crosses indicate death of 
all 5 animals. All Smed-notl (RNAi) animals are able to produce 
blastema cells, independent of the day of transection (B-F, day 4 of 
regeneration). However, the size of the blastema generated 
strongly depends on the day of transection. Animals cut earlier 
produce larger blastemas. Animals cut only 1 day after RNAi are 
able to regenerate photoreceptors (B, day 8 of regeneration) 
although later than control(RNAi) animals (A, day 6 of regenera- 
tion). All blastemas produced by Smed-notl (RNAi) animals 
eventually regress (B-F). (G-H) Intact control(RNAi) (G) and Smed- 
notl (RNAi) (H) animals 20 days after RNAi, anterior side is to the 
top. Smed-notl animals 20 days after RNAi display variable levels of 
head regression defects. Scale bars: 500 um. 
(TIF) 

Figure S4 FACS analysis of planarian cell populations in Smed- 
notl (RNAi) animals. (A- B) FACS profiles of planarian cell 
populations in Smed-notl (RNAi) animals, control(RNAi) animals 10 
and 15 days after RNAi (A) and animals 24 hours after irradiation 
(B). Planarian cells are dissociated and separated by FACS using a 
nuclear dye (Hoechst) and a cytoplasmic dye (Calcein). For RNAi 
animals, two biological replicates were technically replicated twice. 
Similarly, irradiated animals were technically replicated. Gating 
conditions to analyse percentage of XI cells are indicated. Smed- 
notl (RNAi) animals show a mild but significant decrease in 
percentage of XI cells (A, lower row), while irradiation almost 
completely eliminates XI cells (B). 
(TIF) 

Figure S5 Dynamics of neoblasts and their progeny in Smed- 
notl (RNAi) and Smedwi-2(RNAi) animals. (A-F) WMISH of the 
neoblast marker Smedwi-1 (A, D), the early neoblast progeny 
marker Smed-nb.21.1 le (B, E) and the late neoblast progeny marker 
Smed-agat-1 (C, F) in control(RNAi) animals (A-C) and Smed- 
notl (RNAi) animals (D-F) 20 days after RNAi. The level of 
Smedwi-1 signals in Smed-notl (RNAi) animals is variable, including 
animals with almost normal expression (D, top panel) and animals 
with a prominent reduction in Smedwi-1 levels (D, bottom panel). 
All Smed-notl (RNAi) animals present a severely reduced number of 
Smed-nb.2 1 .1 le-positive cells (E). The number of Smed-agat-1- 
positive cells is also variable (F), but all animals have reduced 
levels in the anterior part, typical behaviour of the marker Smed- 
agat-1 upon neoblast perturbation. (G— R) WMISH of the neoblast 
marker Smedwi-1 (G, J, M, P), the early neoblast progeny marker 
Smed-nb.21.1 le (H, K, N, QJ and the late neoblast progeny marker 
Smed-agat-1 (I, L, O, R) in control(RNAi) animals (G-I) and Smedwi- 
2 (RNAi) animals 10 J-L), 15 (M-O) and 20 (P-R) days after 
RNAi. Smedwi-2(RNAi) animals have detectable expression of 
Smedwi-1 in almost all time points (J, M), although a severe decline 
in the level of Smedwi-1 signals is detected 20 days after RNAi (P). 
The dynamics of progeny markers is also abnormal, with a 
progressive decline o£ Smed-nb.21.1 le signals (N, Q) and of Smed- 
agat-1 signals (L, O, R) that precedes the neoblast loss. Anterior is 
to the left. Scale bars: 500 urn. 
(TIF) 

Figure S6 Dynamics of stem cell transcripts and progeny 
transcripts after irradiation. (A-B) Quantification of the level of 
expression by qRT-PCR of the stem cell markers Smedtud-1, 
Smedvas-1, and Smed-pcna (A) and of the neoblast and progeny 
markers Smedwi-1, Smed-nb.21.11e and Smed-agat-1 (B) in animals 1, 
3 and 5 days after irradiation, normalized expression and relative 
to non irradiated samples. Error bars represent standard deviation. 
Animals 1 day after irradiation have around 10% of Smed-pcna 
transcripts of non-irradiated controls, and this number further 
decreases 3 and 5 days after irradiation. However, the expression 
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of Smedtud-1 and Smedvas-1 mRNAs only decreases to around 70% 
and 40% respectively of the level of non irradiated controls, 
reflecting expression that does not localize to neoblasts and is 
therefore not eliminated by irradiation. Similar to Smed-pcna, the 
level of Smedwi-1 transcripts decreases to around 10% of the 
expression in non irradiated controls and becomes almost 
undetectable later. The levels of the progeny specific mRNAs 
Smed-nb.21.11e and Smed-agat-1 decrease progressively at later time 
points of irradiation. Therefore, around 90% of the neoblast 
specific transcripts are eliminated only 1 day after irradiation while 
most of the expression of progeny specific transcripts is still 
detected and the non-neoblast expression of Smedtud-1 and Smedvas- 
1 localized in the CNS is not eliminated by irradiation. 
(TIF) 

Table SI In silico search of CCR4-NOT complex components in 
S. mediterranea. Summary of the CCR4-NOT complex components 
found in silico in S. mediterranea. Each of the described components 
of the yeast, Drosophila melanogaster and human CCR4-NOT 
complexes is indicated. The column "S. mediterranea" indicates 
the given name for each of the genes, the column "Contig Smed 
genome" indicates the genomic contigs in which each locus was 
found, followed by the transcriptomic datasets from Blythe et al. 
and Adamidi et al. The column "D. japonica" indicates the 
accession numbers of the CCR4-NOT components previously 
described in this planarian species. Nine different components of 
the CCR4-NOT complex were found both in genomic and 
transcriptomic sequences, corresponding to the orthologues of all 
metazoan CCR4-NOT complex components. Similar to humans, 
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